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1. Introduction 

In recent years much progress has been made in our understanding of the strongly-correlated 
low energy regime of QCD in terms of gauge-fixed correlation functions. This progress is tightly 
linked to the advance in our understanding of the basic phenomena of low energy QCD, strong 
chiral symmetry breaking and confinement. Confinement denotes the absence of colour charged 
asymptotic states. For infinitely heavy quarks this can be related to the free energy of a single 
quark F q : Bringing a single quark into the system would require an infinite amount of energy in the 
confined phase, whereas a finite energy would suffice in the deconfined phase. The operator which 
creates a static single quark is related to the Polyakov loop L, 

L = —tifPi?*f° d *> A oM , and (L) ~ exp- F " /T , (1.1) 

with the trace in the fundamental representation of SU (N c ) of the closed Wilson line in time di- 
rection. Here, V denotes path ordering, A is the gauge field and the inverse temperature j8 = l/T. 
This links the confinement-deconfinement phase transition in the Yang-Mills system to the order- 
disorder phase transition of center symmetry Zn in SU (N c ): Under center transformations z e Zjv the 
Polyakov loop transforms with L -*■ zL in the fundamental representation. In the center-symmetric, 
confining phase the order parameter vanishes, (L) = 0, while in the center-broken, deconfined phase 
it is finite, (L) * 0. 

It has been shown that L[(Ao)] also constitutes an order parameter for static quark confinement 
[2, 3], as L[(Ao)] > (L[Ao]) in the deconfined but strictly L[(Ao}] = in the confined phase. How- 
ever, L[(Ao)] is easily accessible within functional continuum methods, as (Aq) can be straightfor- 
wardly computed from the equation of motion, that is via the full effective potential V[Ao]. It was 
shown within the functional renormalisation group (FRG) approach [2, 3, 4, 5], that V^Aq] has a 
closed representation in terms of the full propagators of the gluons, ghosts (and quarks) in constant 
Ao-backgrounds. With this connection a confinement criterion in terms of the infrared behaviour 
of the propagators was derived [2]. 

By now (L[Ao]) and the effective Polyakov loop potential V[Ao] has been computed in Yang- 
Mills theory in the Landau gauge [1, 2, 6] and in the Polyakov gauge [3]. The potential obtained 
in different gauges agree well which provides a check of the gauge independence. Recently, these 
studies have been extended to the Coulomb gauge [7] . More recently, lattice results for the potential 
have been obtained [8, 9]. In [4] the Yang-Mills studies have been extended to fully dynamical two- 
flavour QCD. In [10, 11], V[Ao] is used in an effective model approach to QCD with interesting 
applications to thermodynamic observables. 

Here we report on our work in [1], in which we compute the Polyakov loop effective potential 
V[Ao] in the background field formalism in Landau-deWitt gauge within different functional ap- 
proaches, [1]. In Section 2 we derive representations for V[Ao] with the FRG, Dyson-Schwinger 
equations (DSEs) and two-particle irreducible (2PI) functionals. In Section 3 we deduce a confine- 
ment criterion in terms of the propagators: Infrared suppression of gluons but non-suppression of 
ghosts confines static quarks at sufficiently low temperatures. In Section 4 numerical results for 
V[Ao] are presented, which are based on the finite temperature propagators computed within the 
FRG in [12, 13, 14]. 
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d t T k [A;<j>] = 




(a) Flow equation for the Polyakov loop potential. (b) DSE for the background gluon one-point function. 

Figure 1: Functional equations for the (derivative of) the Polyakov loop potential. The gluons are repre- 
sented by wiggly, the ghosts by dashed lines. Filled circles indicate full correlators. Crossed circles denote 
the regulator insertion d t Rk. 



2. Polyakov Loop Potential from Functional Methods 

Most functional approaches are based on closed expressions for the effective action or deriva- 
tives thereof in terms of full correlation functions. Hence, the knowledge of the latter in constant 
Ao-backgrounds allows us to compute the Polyakov loop potential V[An]. In turn, confinement 
requires V[Ao] to have minima at the confining values for Ao, i.e. in SU (N c ) the center-symmetric 
points. With (c)c being the (anti-)ghosts, in the background formalism the Polyakov loop po- 
tential is the gauge-invariant effective action, T[Ao;0] with <p - (a - A-Aq,c,c), evaluted on a 
constant Ao-background and vanishing fluctuation fields, V[A™ nst ] := r[A™ nst ;0]/ (j8V). V is the 
three-dimensional spatial volume. The Polyakov loop, (1.1), is then computed at the minimum 
(Ao) := Ao, m in- Functional equations for the effective action depend on the correlation functions of 
fluctuation fields only. By gauge covariance the correlation functions in the background Landau 
gauge, r(")[A], are directly related to those in the Landau gauge, r(")[0]. Neglecting subleading 
corrections, for more details see [1], the two-point functions of gluon and the ghost on constant 
backgrounds are schematically given by 

d, 2) [A;0] = £ {-D 2 )Z LIT P L I T , and if >[A;0] = {-D 2 )Z C , (2.1) 

L/T 

where the P L I T are the projection operators for the chromoelectric/chromomagnetic modes, and 
Z L / T f c (-D 2 ) ,-D ) are the wave-function renormalisations of the corresponding gluon and ghost 
modes evaluated on covariant momenta D. The FRG representation of the Polyakov loop potential 
has been used in [2, 3,6]. The exact flow equation for the effective action, T^fA; 0], is given solely 
in terms of the propagators, 

d,r*[A^]=jf(i(G a )^ 

where the integration and Matsubara summation involve the respective modes. Further, t = Ink, 
with k being a infrared cut-off scale controlled by the regulator R a i c . At k the propagators are 

G a i c , k {p) = {{T a i,) k {p)+R k {p)y l . (2.3) 
The diagrammatic representation of the flow equation is given in Fig. 1(a). 
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Finally, by integration we get the Polyakov loop potential 

V[A Q ] = V k=A [A ] + -L jf dtd t T k [A ;0] , (2.4) 

where V^' =A [Ao] -> for A/r » 1. The DSE for Yang-Mills theory relevant for the Polyakov loop 
potential is that originating in a derivative of the effective action T w.r.t. the Ao-background at fixed 
fluctuation. It can be written in terms of full propagators and vertices and the renormalised classical 
action, 

8(r[A Q ;0]-S A [A ;0]) _ 1 (3) _ (3) 1 (4) 3 (3) (4) 2 (3) 

with j, and , denoting derivatives of the renormalised classical action 5a and effective 
action with respect to the fields </>,-. For constant backgrounds the DSE equation for Polyakov loop 
potential is given by 

dV[A ] 1 dr[A ;0] 
<9A jSV <9A ' ^ ' ; 

with (2.5) for the right hand side of (2.6). In comparison with the DSE for fluctuating fields only 
the vertices differ. This also leads to the emergence of the two-loop diagram involving ghosts. 
Diagrammatically the DSE is given in Fig. 1(b). 

For details on the different representations see [1]. There it is also shown that for the present 
purpose the 2PI representation of the Polyakov loop potential is identical to the DSE (2.5). 



3. Confinement Criterion 

Here we discuss the restriction of the infrared behaviour of the propagators that follows from 
the FRG. Here, all loops in the full flow, (2.2), effectively boil down to expressions of the form 

. M yfii^ wo (31) 

,^7o (2?r)3 xZ{-Dlp 2 ) + ^{DM)+Rk{x)' 

where x = -Dq + p 2 and Z = Zi,Zj,Z c . In (3.1) we do not include the overall minus sign of the 
ghost contribution. The function stands for the scalar parts of the regulators, which are chosen 
to be Z-independent and 0(4)-symmetric, i.e. Rk(x) -xr(x/k 2 ). The terms Am 2 (D,0) = 0, with 
Am 2 = Am 2 , Am 2 , Am 2 , are thermal corrections to the Polyakov loop, which disappear with powers 
of the temperature for T -*■ 0, and for the present purpose we will neglect these terms, [1], 

Am 2 (D,A Q ) = 0. (3.2) 

At vanishing temperature we regain 0(4)-symmetry and Z{-D^,p 2 ) -> Z(-Dq + p ). In terms of 
scaling coefficients K at vanishing temperature, 

Z(x)=Zx K , (3.3) 
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the corresponding propagators exhibit non-thermal mass gaps m gap for K < 0. Computing the ex- 
pression (3.1) leads to a deconfining potential in Aq: Taking a derivative w.r.t. A$ it vanishes at 
Ao = and the center-symmetric points, i.e. where the Polyakov loop vanishes, L = 0. The second 
derivative is positive at Ao = while it is negative at the center-symmetric points. 

Hence, with (3.2) the gluon loops in the flow equation (2.2) give deconfining potentials for 
the gluons. In contrast, the ghosts give a confining contribution due to the relative minus sign. 
Thus, confinement in covariant gauges requires the suppression of gluonic contributions. Indeed, 
the same mechanism is at work in the DSE representation [1]. With the assumption (3.2) this yields 
a necessary condition for confinement, 

lim — — = and lim — — > , thus, K L/T < and K c > , (3.4) 

*-0Z L/r (x) x-*0Z c (x) 

where K L i T are the anomalous dimensions of the chromoelectric and chromomagnetic gluon. In 
(3.4) the condition involving the ghost guarantees that its contribution dominates the trivial gauge 
mode. Eq. (3.4) is a sufficient condition for confinement in Yang-Mills propagators (with the 
assumption (3.2)) as it leads to a vanishing order parameter, the Polyakov loop expectation value. 
Note also that it has been shown [15, 16] within a scaling analysis that > has to hold in the 
infrared. 

The infrared limits (3.4) are satisfied for both decoupling type, i.e. K c = and K a - - I, and 
scaling type, i.e. k c -K - - K a /2 with K w 0.595, solutions. For an overview of the different solutions 
see e.g. [17]. As a consequence, both types exhibit confinement. Furthermore, they do not affect 
the critical physics at the phase transition, as it is triggered in the mid-momentum regime where 
both types of solutions agree on a quantitative level. 

The confinement criterion provided here a leads to simple counting scheme for general theories 
and gauges, e.g. Yang-Mills-Higgs systems [18, 19, 20, 21], or QCD with fundamental or adjoint 
quarks in the Landau gauge, [4, 5, 22, 23]. The respective scenarios are detailed in [1]. 



4. Results for the Polyakov Loop Potential and T c 

Here we present the results for V[Ao] in the FRG, DSE and 2PI approaches, which have 
previously been computed with the thermal propagators obtained from the FRG formalism [12, 13, 
14]. The potential determined from the flow equation (2.2) is depicted for temperatures around the 
confinement-deconfinement phase transition for SU(2) in Fig. 2(a), and for SU(3) in Fig. 2(b). 
The critical temperatures T c are determined at 

T c su( - 2) =230±23MeV and T c su( - 3) = 275 ±27MeV. (4.1) 

The absolute temperatures in (4. 1) are set in comparison to absolute scales of the lattice propagators 
[24, 25, 26] with a string tension a = (420 MeV) 2 , which leads to the dimensionless ratios of 
Tc U ^ 2 l\/<5 ~ 0.548 and Tc U /y/o & 0.655, respectively. For SU (3) our result is in quantitative 
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agreement with the lattice results, see e.g. [27], 

j,SU{2) T SU(3) 

T su(2) = 295 MeV T su(3) = 270Me v thus, w 0.709 , w 0.646 . (4.2) 

t,iatt c.iatt 

Our relative normalisation is taken from the peak position of the lattice propagators in [24, 25, 26]. 
This normalisation has an error of approximately 10%, which sets the errors. 

For SU (2) the critical temperature r c 5C/ ( 2 ) is significantly lower than the lattice temperature in 
(4.2). This is linked to the missing back-reaction of the fluctuations of V[Ao] onto the propagators 

SU(2) 

[1]. Its inclusion in the Landau gauge flow leads to a transition temperature of T c = 300 ± 
30 MeV [28], in quantitative agreement with results in Polyakov gauge [3] and from lattice gauge 
theory. 

For the computation of V[Ao] from its DSE, (2.6), two-loop diagrams involving dressed ver- 
tices in the presence of a background Aq need to be determined. However, employing a resumma- 
tion scheme based on the 2PI-hierachy their contributions can be minimised, which sets the optimal 
renormalisation point jj. w lGeV [1]. As a consequence, it suffices to evaluate a one-loop equation 
which does not involve dressed Green functions other than the propagators. 

The one-loop diagrams in the DSE, Fig. 1(b), reduce to 

V'(q>)= l -V^M + flKT{n + q>) 

x[ Za G L (q + 2nT(p) + 2z a G T (q + 27lT(p) -2z c G c (q + 27tT<p)] , (4.3) 

where G L / T , G c are the propagators of the chromoelectric/chromomagnetic gluon and ghost nor- 
malised at z a / c = Z a / c (p 2 = pL 2 ) at zero temperature. The first term originates in the trivial gauge 
mode. The evaluation of (4.3) leads to the potential depicted in Fig. 3(a) for SU{2) and Fig. 3(b) 
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(a) Potential for SU(2). (b) Potential for 5(7(3). 

Figure 3: Polyakov loop potential for SU(2) and SU(3) from its DSE/2PI-representation within an opti- 
mised RG-scheme. 



for SU (3). From those the confinement-deconfinement critical temperatures are determined at 



,su(2) T su(3) 



T c su{2) = 235±24MeV, T c su{3) = 275 ±27MeV, thus, 




* 0.655. (4.4) 



5. Conclusions 

We have studied the confinement-deconfinement phase transition of static quarks with the 
Polyakov loop (effective) potential, which was computed in the background field formalism in the 
Landau-deWitt gauge [1]. This allowed for the use of Landau gauge propagators at vanishing 
background, which had been obtained previously in the framework of the FRG [13, 12, 14]. The 
study here was done within different non-perturbative functional continuum approaches, i.e. using 
FRG-, DSE- and 2PI-representations of the effective action. 

In agreement with lattice gauge theory we find a second and first order phase transition for 
SU (2) and SU (3), respectively. The critical temperatures in different methods agree quantitatively 
within errors. These results are stable w.r.t. the choice of scaling or decoupling type solutions for 
the propagators. The critical temperature for SU (3) from functional methods agrees quantitatively 
with lattice results. For SU (2) the critical temperature from functional methods is below the lattice 
temperature, which is due to the missing back-reaction of the fluctuations of the Polyakov loop 
potential onto the higher Green functions, in particular the chromoelectric propagator. 

The dependence of the Polyakov loop potential on the individual gluon and ghost modes was 
exploited to derive a simple criterion for static quark confinement. For infrared suppressed gluon 
but not suppressed ghost propagators confinement is immanent at sufficiently low temperatures, 
as the transversal gluons do not contribute and the ghost modes overcompensate the deconfining 
longitudinal gluon mode. 

In summary this work establishes the functional approach to confinement with means of the 
Polyakov loop potential. It only requires the background gauge propagators. Hence, it gives easy 
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access to structural questions about the confinement mechanism and to explicit computations of the 
phase structure. 
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